Hall-Lorenz number paradox in cuprate superconductors 
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Significantly different normal state Lorenz numbers have been found in two independent direct 
measurements based on the Righi-Leduc effect, one about 6 times smaller and the other one about 2 
times larger than the Sommerfeld value in single cuprate crystals of the same chemical composition. 
The controversy is resolved in the model where charge carriers are mobile lattice bipolarons and 
thermally activated nondegenerate polarons. The model numerically fits several longitudinal and 
transverse kinetic coefficients providing a unique explanation of a sharp maximum in the temperature 
dependence of the normal state Hall number in underdoped cuprates. 

PACS numbers: 74.25.Fy,74.20.-z,72.15.Eb,74.72.-h 



Particular interest in studies of high-temperature 
superconductors lies in a possible violation of the 
Wiedemann- Franz (WF) law in doped cuprates. A de- 
parture from the Fermi/BCS liquid picture is seen in both 
the superconducting and normal state thermal conduc- 
tivities and might be related by a common mechanism 1, 
2\. Takenaka et al.Q systematically studied the oxygen- 
content dependence of the insulating state thermal con- 
ductivity enabling them to estimate the phononic con- 
tribution, K p h(T) for the metallic state to some extent. 
Their analysis led to the conclusion that the electronic 
term, k is only weakly T-dependent. This approximately 
T-independent k in the underdoped region therefore im- 
plies the violation of the WF law since the resistivity 
is found to be a non-linear function of temperature in 
this regime. A breakdown of the WF law has been seen 
in other cuprates such as Pr2~ x Ce x CuO\ at very low 
temperatures]!! . On the other hand measurements by 
Proust et al.|3] on T12Bo,2CuOq + s have suggested that 
the Wiedemann-Franz law holds perfectly well in the 
over doped region. However in any case the extraction 
of the electronic thermal conductivity has proven diffi- 
cult and inconclusive as n and n p h are comparable at 
elevated temperatures, or there is a thermal decoupling 
of phonons and electrons at ultra-low temperatures 4] . 

This uncertainty has been avoided in measurements of 
the Righi-Leduc effect. The effect describes transverse 
heat flow resulting from a perpendicular temperature 
gradient in an external magnetic field, which is a ther- 
mal analog of the Hall effect. Using the effect the "Hall- 
Lorenz" electronic number, Lh = (e/kFj) 2 k xv / (T(7 xy ) 
has been directly measured^ in Y Ba2CuzO§.g$ and 
Y BCL2CU3O6.6 since transverse thermal K xy and electri- 
cal a xy conductivities involve presumably only electrons. 
The experimental Ljy(T) showed a quasi-linear temper- 
ature dependence above the resistive T c , which strongly 
violates the WF law. Remarkably, the measured value 
of Lh just above T c turned out precisely the same as 
predicted by the bipolaron theory 6], L = 0.15Lo, where 
Lq = 7r 2 /3 is the conventional Sommerfeld value. The 
breakdown of the WF law revealed in the Righi-Leduc 
effect p| has been explained by a temperature-dependent 
contribution of thermally excited single polarons to the 



transverse magneto-transport 0. 

Surprisingly more recent measurements of the Hall- 
Lorenz number in single crystals of optimally doped 
YBCI2CU3OQ.95 and optimally doped and underdoped 
EuBa2Cu^O y led to an opposite conclusion The ex- 
perimental Lh for these samples has turned out only 
weakly temperature dependent and exceeding the Som- 
merfeld value by more than 2 times in the whole temper- 
ature range from T c up to the room temperature. Fol- 
lowing an earlier claim[!j Matusiak and WolfQ have ar- 
gued that a possible reason for such significant difference 
might be that Zhang et al.[f| used different samples, one 
for n xy and another one for a xy measurements, which 
makes their results for Lh inconsistent. 

Here I argue that there is no inconsistency in both 
Lh determinations. One order of magnitude difference 
in two independent direct measurements of the normal- 
state Hall-Lorenz number is consistently explained by the 
bipolaron theory ^(j- The theory explains the huge dif- 
ference in the Hall-Lorenz numbers by taking into ac- 
count the difference between the in-plane resistivity of 
detwinned[j| and twinned0 single crystals. The the- 
ory fits well the observed Lh{T)s and explains a sharp 
Hall-number maximum[3 observed in the normal state 
of underdoped cuprates. 

In the presence of the electric field E, the temperature 
gradient VT and a weak magnetic field B|| z _L E and 
VT, the electrical currents in x, y directions are given by 

j x =a xx V x (n ~ 2e0) + a xy V y (fi - 2ecj)) 
j y =ayy\7y(n - 2e<j>) + a yx V x (n - 2ecj)) 

~\~byy^7 yT ~\- hyx^ X^ 1 



(1) 



and the thermal currents are: 

w x =c xx V x (n - 2e(j>) + c xy V y (fi - 2e4>) 

~^~d XX ^/ X T ~\- d X y^ yT 
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CyyVy(fl ~ 2e(f)) + Cy X \7 ' x (jl ~ 2e 



~^~^yy^ ' y^ d~ d yx ^ X T. 



(2) 
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Here /i and are the chemical and electric potentials. 

Real phonons and (bi)polarons are well decoupled 
in the strong-coupling; regime of the electron-phonon 
interactionJlOj so the standard Boltzmann equation for 
the kinetics of renormalised carriers is applied. If we 
make use of the t(E) — approximation [lT| the kinetic co- 
efficients of bipolarons are found as0 (fcs = h = c = 1) 
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where 



(Q(E)) = 



/ °° dEQ(E)ED b (E)df b /dE 
J °° dEED b {E)df b /dE ' 



(3) 



D b (E) oc E 4 / 2 ^ 1 is the density of states of a d- 
dimcnsional bipolaron spectrum, E = K 2 /(2m b ), g b — 
2e/m b , and fb(E) is the equilibrium distribution func- 
tion. Polaronic coefficients are obtained by replacing su- 
per/subscripts b for p, double elementary charge 2e for 
e and /i for /i/2 in all kinetic coefficients, and m b for 
2m p in <Zjj and . The kinetic energy of bipolarons, E 
should be replaced by E + T*, where E = k 2 /(2m p ) is 
the polaron kinetic energy, and T* is half of the bipo- 
laron binding energy (i.e. the pseudogap temperature in 
the theory[n|). 

The in-plane resistivity, p, the Hall number, Rh, and 
the Hall-Lorenz number, Lh are expressed in terms 
of the kinetic coefficients as = 2ea xx , Rh = 

/2eB(a xx ) 2 , and 



C[(dy X Cl XX Cy X b XX ^Q XX C XX (b XX Qy X by X d XX ^)] 
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y x xx 



(4) 



respectively, where a, b,c,d — a p + a b ,b p + b b . c p + c b ,d p + 
d". 

The zn-plane resistivity, the temperature-dependent 
paramagnetic susceptibility, and the Hall ratio have al- 
ready been described by the bipolaron model ta king 
into account thermally activated single polarons^^ Il3t 



E3. The bipolaron model has also offered a sim- 
ple explanation of c-axis transport and the anisotropy 
of cuprates[T3. flfi! IT^ . IT^ |. The crucial point is that 
single polarons dominate in c-axis transport at finite 
temperatures because they are much lighter than bipo- 
larons in c-direction. Bipolarons can propagate across 
the planes only via a simultaneous two-particle tun- 
nelling, which is much less probable than a single po- 
laron tunnelling. However, along the planes polarons and 
inter-site bipolarons propagate with comparable effective 
masses jlO|. Hence in the mixture of nondegenerate quasi- 
two-dimensional (2D) bosons and thermally excited 3D 
fermions, only fermions contribute to c -axis transport, 
if the temperature is not very low, which leads to the 
thermally activated c -axis transport and to the huge 
anisotropy of cuprates ^ • 

We have also shown 3 that by the necessary inclusion 
of thermally activated polarons, the model, Eq.(4) pre- 
dicts a breakdown of the WF law with the small near- 
linear in temperature Hall-Lorenz number, as observed 
experimentally by Zhang et al. 5] (see Fig.l). Let us now 
show that the bipolaron model describes the contrasting 
observations of Ref.|S| as well, if the ratio of bipolaron 
and polaron mobilities, a — 2T b m p /T p m b becomes rela- 
tively small. 

Both polaronic and bipolaronic carriers are not degen- 
erate above T c , so the classical distribution functions, 
f b = ycxp(-E/T) and f p = y^ 2 exp[-(E + T*)/T] are 
applied with y = cxp(/i/T). The chemical potential is 
evaluated using 2n b +n p = x/v^, where x is the number of 
itinerant holes in the unit cell volume vq not localised by 
disorder. The bipolaron density remains large compared 
with the polaron density in a wide temperature range, so 
that n b vo k, xj2 and y « irx / (m b a 2 T) for quasi-2D bipo- 
larons. Then the atomic density of 3D polarons is n p Vo = 
Tm p a 2 exp(-T* /T)(xm p /2ir 2 m b ) 1 / 2 (a is the lattice 
constant). The ratio (3 = n p /2n b remains small at any 
pseudogap temperature T* and any relevant doping level 
x > 0.05, (3 w Texp(~T* /T)(18m p /n 2 xm b )^ 2 /W < 1, 
if the temperature T is small compared with the polaron 
bandwidth W = 6/m p a 2 . Hence, if the mobility ratio a 
is of the order of unity, both longitudinal and transverse 
in-plane magnetotransport is dominated by bipolarons, 
which explains a remarkably low Lu in high-quality de- 
twinned crystals used in Ref. 5], Fig.l. 

On the other hand, twinned crystals used in Ref.0 had 
the in-plane resistivity several times larger than those of 
Ref.0 presumably resulting from twin boundaries and 
long term aging. The twin boundaries and other defects 
are strong scatterers for slow 2D bipolarons (see below), 
while lighter quasi-3D polarons are mainly scattered by 
real optical phonons, which are similar in all crystals. 
Hence one can expect that a becomes small in twinned 
crystals of Ref.||. If the condition a 2 <C (3 is met, then 
only polarons contribute to the transverse electric and 
thermal magnetotransport. It explains about the same 
thermal Hall conductivities ( 10 3 B W/Km 

at T=100K) dominated by polarons in both crystals of 
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FIG. 1: The Hall-Lorcnz number Lh in underdoped twinned 
EuBa2Cv.3O6.65 (circles)|J compared with the theory, Eq.(7) 
when a <C 1 (upper line), and the significantly different Hall- 
Lorenz number in detwinned Y Ba.2Cu3O6.95 (triangles) [fj 
described by the same theoryQ with a moderate value of 
a = 0.44 (lower line). 



Fi3a2Cw.3O6.95 used in Ref.[f| and in Ref.0, and at the 
same time a substantial difference of their electrical Hall 
conductivities, <J xy , as bipolarons virtually do not con- 
tribute to a xy in the twinned samples. 

To arrive at simple analytical results and illustrate 
their quantitative agreement with the experiment @ let 
us assume that a 2 <§; (3, but a > fi, and neglect an en- 
ergy dependence of the transport relaxation rates of all 
carriers. In such conditions bipolarons do not contribute 
to transverse heat and electric flows, but determine the 
in-plane conductivity. Kinetic responses are grossly sim- 
plified as 



Rh = 



v (3 
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L H = 4.75 + 3T*/T + (T*/T) 2 



(5) 

(6) 
(7) 

As in the case of a 2 > /3, discussed in Ref. , the re- 
combination of a pair of polarons into bipolaronic bound 
states at the cold end of the sample results in the break- 
down of the WF law, as described by two temperature- 
dependent terms in Eq.(7). The breakdown is reminis- 
cent of the one in conventional semiconductors caused 
by the recombination of electron-hole pairs at the cold 
end • However, the temperature dependence and the 
value of Lh(T) turn out remarkably different. When 
a 2 <C (3, The Hall-Lorenz number is more than by an or- 
der of magnitude larger than in the opposite regime, a 2 > 
f3. It increases with temperature lowering rather than 



decreases fitting well the experimental observation Q in 
twinned underdoped single crystals of EuBa2Cu^O(,.^ 
with T* = 100K, Fig.l. Hence by varying the bipolaron 
to polaron mobility ratio, a, the model accounts for quali- 
tatively different behaviours of Lh(T) in twinned and de- 
twinned cuprates. The energy dependence of relaxation 
rates might somewhat change numerical coefficients in 
Eq.(7), but it does not qualitatively change the temper- 
ature dependence and the value of Lh{T). 

The Hall-Lorenz number is the ratio of different ki- 
netic coefficients rather than a proper kinetic response 
function. However, its significant departure from the 
Sommerfeld value Lq rs 3.3 clearly indicates a non-Fermi 
liquid behaviour since the relaxation mechanism virtu- 
ally cancels in the ratio. The partially gapped Fermi- 
liquid model used to explain large Lh in Ref. 8f predicts 
a quadratic decrease of Lh(T) with temperature lower- 
ing, rather than a steep increase as observed, Fig.l. To 
account for an unexpected rise of Lh(T) below T « 160K 
in underdoped samples, Matusiak et al.@ suggested an 
opening of another narrower gap. However the gapped 
Fermi liquid model is clearly incompatible with the near 
temperature-independent resistivity and with the sharp 
maximum of the normal state Hall ratio at 100K, as mea- 
sured in Ref. 8], Fig. 2. It is also hard to accept the claim 
of Refs. [a, Oj that the research team of Ref. |£j could so 
badly manipulate their data to arrive at an erroneous 
Lh more than one order of magnitude smaller in identi- 
cal cuprates. 

On the contrary our model explains the near 
temperature- independent resistivity and the unusual Hall 
ratio, Fig. 2. If we assume that in EuBa,2Cu? ) 0§.§§ slow 
bipolarons are mainly scattered by neutral defects and 
twin boundaries, their relaxation rate depends on the 
temperature as Tbo/ri, = l + (T/To) 1 / 2 , where no is a con- 
stant. The temperature independent contribution comes 
from the scattering rate off neutral impurities with the 
carrier exchange 11] similar to the scattering of slow elec- 
trons by hydrogen atoms. The square-root term origi- 
nates in the scattering of slow bipolarons by point de- 
fects and twin boundaries with a temperature indepen- 
dent mean-free pass. The scale T thus depends on the 
relative strength of two scattering mechanisms. The the- 
oretical resistivity 



P(T) 
Po 



1 + (T/Tq) 1 / 2 



(8) 



fits well the experimental p(T) in the entire normal- 
state region with po = rn^VQ / \2e 2 xt(,o) — 1.3x 
10" 5 fim and T =321 K, Fig.2. Lighter 3D polarons 
are scattered by defects and optical phonons, so that 
Tpo/Vp = {T/W) 1 / 2 + £?exp(— uj/T) with a temperature- 
independent T p o- Then, using n p oc Texp(—T*/T), 
Eq.(6) yields 



R H (T)=p 2 (T)- 



ATexp(-T*/T) 



(9) 



[T 1 /2 + 6 CX p(_ w /r)]2- 

This expression fits extremely well the experimental 
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FIG. 2: The in-plane resistivity p (triangles) and the Hall 
ratio Rh (circles) of underdoped twinned EuBa2Cus.OQ.Qz,, 
Ref.0, compared with the theory (lines). 



Rh(T) with temperature independent constants A = 
e 3 T 2 (18xm p /T: 2 m b ) 1 / 2 /(vom 2 )=275 m/tt 2 C and b = 

BW 1 / 2 =122K 1 / 2 , the reasonable value of the character- 
istic optical phonon frequency u — 470K, and the same 
pseudogap T* = 100K as in the Hall-Lorenz number in 
Fig.l. It appears almost perfect even in the critical re- 
gion very close to T c , Fig. 2, if one uses the experimental 
p(T) in Eq.(9). However the maximum of the Hall ratio 
is a normal state feature lying well above the critical re- 



gion by about 30K, Fig. 2, as in underdoped Y Ba2Cu^O y , 
Ref . jlflj . At temperatures below T* the Hall ratio drops 
as the number of thermally activated polarons decreases, 
and at temperatures above T* it drops since the polaron 
relaxation time decreases. 

To verify the self-consistency of the model let us esti- 
mate a and /?. In the optimally doped samples one ex- 
pects a 2 of the order of (3, so the Hall ratio approximately 
measures the itinerant carrier density, i?^ 4 ks v /ex opt - 
Then using the experimental values[2j of Rjj in optimally 
doped EuBa2CuzOi and underdoped EuBaiCu-$0§.§§ 
one estimates a 2 / ~ R'S 'x op t / 'xRh J$ 0.1 in under- 
doped EuBa.2Cu3O6.65, which justifies one of our as- 
sumptions. To get a > (3 we have to assume that (5 < 0.1, 
which is indeed the case in the whole temperature range, 
if the polaron band is wide enough, W > 5000K. Fi- 
nally using the values of po and A and taking x = 0.1, 
nip = 5m e , vq = 0.2 nm 3 and ra b = 2m p the polaron and 
bipolaron mean-free pass is estimated as l p « 4 nm and 
lb ~ 0.3(m(,/m e ) 1 / 2 nm, respectively (here m e is the free 
electron mass) . Their values are large compared with the 
lattice constant justifying the Boltzmann approximation 
for all carriers. 

To sum up, the bipolaron theory resolves the paradox 
of very different Hall-Lorenz numbers found in two inde- 
pendent measurements ja, la] in cuprate single crystals. It 
explains a flat temperature dependence of the in-plane 
resistivity and the sharp maximum in the normal-state 
Hall number of underdoped cuprates as well. 
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SRC (UK) (grant EP/C518365/1) and enlightening dis- 
cussions with Nigel Hussey of thermal conductivity mea- 
surements. 
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